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As undesirable defects, grain boundaries (GBs) are widespread in epitaxial graphene using existing 
growth methods on metal substrates. Employing density functional theory calculations, we first 
identify that the misorientations of carbon islands nucleated on a Cu(lll) surface lead to the 
formation of GBs as the islands coalesce. We then propose a two-step kinetic pathway to effectively 
suppress the formation of GBs. In the first step, large aromatic hydrocarbon molecules are deposited 
onto a v3 x y/3 superstructured Cu-Mn alloyed surface to seed the initial carbon clusters of a 
single orientation; in the second step, the seeded islands are enlarged through normal chemical 
vapor deposition of methane to form a complete graphene sheet. The present approach promises to 
overcome a standing obstacle in large scale single-crystal graphene fabrication. 



PACS numbers: 81.05.ue, 81.15.-z, 68.43.-h, 68.35.bd 

Graphene is a one-atom-thick flat sheet of carbon 
atoms packed into a honeycomb structure. Because of 
its superb mechanical, electronic, optical, and thermal 
properties [I] [2J , graphene has limitless potential for fu- 
ture device applications. To fully realize the functionality 
of graphene, it is highly desirable to fabricate large-scale 
monolayer graphene with no or minimal structural de- 
fects. Among the different fabrication avenues being ex- 
plored, epitaxial growth on transition metal substrates 
using hydrocarbon or other carbon sources stands out 
as a highly appealing approach [5HTT]. especially on Cu 
surfaces [IT] . Cu has the merit of low carbon solubil- 
ity, which leads to a self-limiting growth process con- 
fined to its surface |12j . and diverse carbon sources can 
be used to grow graphene on Cu substrates [T51 - fT5] . The 
relatively weak carbon-copper interaction compared to 
carbon-carbon interaction enables fast diffusion of carbon 
atoms and efficient nucleation of carbon islands across 
the whole surface [IB], indicating the feasibility of mass 
production of epitaxial graphene. Indeed, it has been re- 
ported recently that the single-crystal domains of mono- 
layer graphene grown on Cu can reach the dimensions of 
0.5mm on a side |17j . Furthermore, transfer of epitaxial 
graphene to other substrates can be readily achieved via 
chemical etching jTBJ [T9] . 

However, one standing obstacle facing the commu- 
nity of epitaxial graphene on Cu is the prevalence of 
grain boundaries (GBs) undesirably introduced during 
growth [T9- 26J. A grain boundary refers to the junc- 
tion region of two crystalline grains with different ori- 
entations. The detailed atomic structures of the GBs 
in epitaxial graphene have been investigated exten- 
sively [T9J l2ITf2o] , and their presence has been shown to 



severely degrade the electronic, transport, and mechani- 
cal properties of graphene [19j [24] • Experimental efforts 
have also been made to suppress their creation during 
growth [T7J [T9J [27] , but so far with limited success, partly 
because the underlying formation mechanism of the GBs 
is still unclear. Existing experimental observations sug- 
gest that the GBs can form in the initial nucleation stage 
when several graphene grains emanate from one nucle- 
ation site [23]; alternatively, they can be formed in the 
later growth stage when different graphene grains with 
relative misorientations coalesce [T3J [HI US] ■ 

In this letter, we first demonstrate that, because of 
the inherently weak C-Cu interaction, orientational dis- 
orders of carbon islands on Cu(lll) will be abundant in 
the early stages of nucleation and growth. Such disor- 
ders cannot heal themselves with the enlargement of the 
islands, leading to the prevalence of graphene GBs upon 
island coalescence. Based on this understanding, we pro- 
pose to use a functionalized Cu(lll) surface to lift the 
energy degeneracy in the early stages of nucleation and 
growth, thereby suppressing orientational disorders of the 
islands and the subsequent GBs. Our proposed kinetic 
pathway invokes the steps of "seed and grow" [25J . In the 
seeding step, carbon clusters are initiated by depositing 
coronene [29] on a (-y/3 x -y/3) R30° Mn-Cu(lll) alloyed 
surface (SUl EI] , which effectively helps the islands to se- 
lect predominantly only one orientation on the super- 
structurally alloyed surface. In the growing step, larger, 
monolayer graphene is formed by conventional chemical 
vapor deposition (CVD). 

Our density functional theory (DFT) calculations are 
carried out using the Vienna ab initio simulation package 
(VASP) [32] with PAW potentials [Ml 133 and the gen- 
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FIG. 1: (color online), (a) and (b) Structural illustrations of 
the Cu (111) substrate and graphene, where the yellow and 
red dashed lines show their respective high-symmetry axes, 
(c) and (d) Illustrations of two geometries where a 7CR car- 
bon cluster is at a HSO on the Cu(lll) surface. In (c), the 
edge C atoms reside at the 3-fold hollow sites; in (d), the edge 
C atoms are at the bridge sites between two surface Cu atoms. 

eralized gradient approximation (PBE-GGA) [35] for the 
exchange-correlation functional. The lattice constant of 
Cu is obtained via structural optimization. The generic 
Cu(lll) surface is modeled by a slab of 5 atomic lay- 
ers, and the Mn-Cu(lll) surface is realized by substitut- 
ing Cu atoms by Mn atoms at ordered positions in the 
first layer. The vacuum layers are more than 13 A thick 
to ensure decoupling between neighboring slabs. Dur- 
ing relaxation, atoms in the lower 2 atomic layers are 
fixed in their respective bulk positions, and all the other 
atoms are allowed to relax until the forces on them are 
smaller than 0.01 eV/A. A 2x2x1 fc-point mesh is used 
for the 6x6 surface unit cell and 3x3x1 for the 4x4 sur- 
face unit cell [3B]- The calculations with Mn atoms are 
spin-polarized. We consider the ferromagnetic configu- 
ration of the Mn-Cu(lll) surface, because of the trian- 
gular arrangement of the Mn atoms. The binding en- 
ergies are calculated as AE = S adsorba t e + £ substrate - 

-^adsorbatc+substrate • 

Crystalline Cu has a face-centered cubic (fee) struc- 
ture, and its (111) surface exhibits a hexagonal packing 
of surface atoms. As shown in Figs. 1 (a) and (b), both 
the Cu(lll) surface and graphene have atomic arrange- 
ments with six-fold symmetry. Therefore, if carbon clus- 
ters nucleated at different sites are all oriented at the 
same high-symmetry orientation (HSO) of the Cu(lll) 
surface (e.g., Fig. 1 (c) or (d)), their structural coher- 
ence will be ensured by the Cu substrate and there will 
be no GBs when they merge. However, when a sim- 
ple six-fold symmetric carbon cluster composed of seven 
6-carbon rings (7CR) is placed on the Cu(lll) surface, 
our detailed DFT calculations reveal that the energet- 




FIG. 2: (color online). Side and top view of a 7CR on the 
Cu(lll) surface, illustrating the domed nature (a) and the 
rotated nature (b) from the HSO of Fig. 1 (c), respectively. 

ically most stable geometry deviates from the HSO of 
Fig. 1 (c) by 11° (Fig. 2 (b)). The carbon cluster also 
has a dome-like structure (Fig. 2 (a)), with the central 
C atoms ~2.30 A from the Cu surface. Therefore, the 
cluster remains strongly bonded to the substrate only at 
the periphery while the interaction between the central 
C atoms and the substrate is rather weak, similar to the 
domed structure on Ir(lll) [37]. Each of the 12 edge 
C atoms has two C neighbors, and prefers to reside at 
the bridge sites between two surface Cu atoms, because 
these edge atoms are closer to sp 3 hybridization than sp 2 , 
thereby providing the driving force for the rotation of the 
island away from the HSO. 

Now we go from the early stages of island nucleation 
and growth to island enlargement and coalescence. Since 
a small cluster such as the one shown in Fig. 2 is not 
oriented at the HSO, there will be a degenerate mirror 
geometry with respect to the symmetry axis of Cu(lll), 
indicating that islands with relative misorientations can 
coexist. As a cluster grows larger, more edge C atoms 
will be involved in determining its preferred orientation 
by adjusting their bonds with the underlying Cu atoms. 
Therefore, there will be more nearly degenerate orienta- 
tions, thus broadening the orientational disorder of the 
carbon clusters. When a cluster has grown large enough 
such that the edge atoms contribute only minimally to 
the total binding energy, the cluster is either still in an 
energetically stable orientation different from the HSO, 
or is too large to adjust its orientation to an energetically 
more favorable HSO. When two such clusters with a rel- 
ative misorientation coalesce, a larger island containing a 
GB is formed, with a characteristic angle defined by the 
initial misorientations of the merging clusters and the 
local structural adjustment within the boundary. This 
scenario is qualitatively consistent with existing experi- 
mental observations, and the detailed distribution of the 
GB angles may also depend on the specific growth con- 
ditions |2"5H2"5] . 

Next we search for ways to lift the energy degener- 
acy in island nucleation and growth via Cu(lll) surface 
modification. Fig. 3 (a) illustrates a (y3 x v3) R30° 
X-Cu(lll) superstructurally alloyed surface, which has 
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FIG. 3: (color online), (a) Structure of the (V3 x \/3) R30° 
Mn-Cu(lll) superstructurally alloyed surface, (b-d) Three 
stable geometries of a 7CR adsorbed on the Mn-Cu(lll) sur- 
face. As indicated by the red dashed lines, (b) and (c) il- 
lustrate two HSO configurations, while (d) is rotated from a 
HSO, and their relative stabilities are indicated by their total 
energy differences. 

transition metal atoms X substituting Cu at ordered po- 
sitions. The choice of X is guided by the requirements 
that, (a) C binds more strongly to X than to Cu, so that 
the nucleated carbon clusters will prefer a HSO in max- 
imizing their interaction with the X atoms; and (b) the 
alloyed X-X atoms are repulsive. We find that a num- 
ber of transition metals (X= Ru, Fe, Co, Ni, Mn) have 
unfilled d orbitals and bind more strongly with C, but 
only Mn atoms are mutually repulsive when alloyed into 
the top layer of Cu(lll) [35], indicating that Mn could 
form a superstructured surface alloy with Cu. Indeed, ex- 
perimentally the formation of a (y/3 x \/3) R30° super- 
structured Mn-Cu(lll) surface has been observed [3TI] . 
whereas the other metals only form islands or overlayers 
on Cu surfaces [39] . 

Still choosing 7CR as the testing baby graphene, we 
then calculate the total energies of 7CR with different 
orientations on the (y/3 x y/S) R30° Mn-Cu(lll) alloyed 
surface. We find three stable or metastable configura- 
tions of a 7CR island through structural optimization, 
differentiated by placing the center of the 7CR above 
a Mn or Cu atom, as shown in Fig. 3. Two of them 
(Figs. 3 (b) and (c)) are at HSO, but only the HSO in 
Fig. 3 (b) is the most stable, while the energy of the 
other two configurations is higher by 0.32 eV and 0.45 
eV, respectively. Therefore, the Mn atoms alloyed into 
the Cu(lll) surface indeed successfully help to pin the 
7CR at the HSO. To see the underlying atomistic rea- 
son, we note that in all the three cases, the island has 
a dome-like geometry similar to that on a pure Cu(lll) 



surface, indicating the predominant interaction with the 
substrate at its edge. Moreover, in the most stable con- 
figuration, the 7CR maximizes its contact with the Mn 
atoms at the periphery. The calculated binding energy 
per edge C atom of 7CR is 0.63 eV on Cu(lll) and 0.89 
eV on Mn-Cu(lll). 

To take advantage of the superstructural Mn-Cu(lll) 
surface and effectively suppress the possible disorders in- 
duced in the initial nucleation process, we propose the 
use of coronene as a good candidate to seed the initial 
carbon clusters. As a polycyclic aromatic hydrocarbon, 
coronene [29] is just like the 7CR island, but with a hy- 
drogen atom on each edge C atom. We have compared 
the dehydrogenation process of coronene to that of ben- 
zene, which has been used for low temperature graphene 
growth on Cu [13]. In free space, the C-H bond dissocia- 
tion energies of benzene and coronene are very close |40j , 
as verified also by our DFT calculations. Considering the 
effect of the catalytic substrates, the energy difference be- 
tween the initial state, where a coronene is adsorbed onto 
the Mn-Cu alloyed surface, and the final state, where all 
the edge hydrogen atoms are dissociated in the form of 
H2, is 1.23 eV per H atom, which is 0.16 eV lower than 
that of benzene on Cu(lll). Therefore, the alloyed sur- 
face will be more catalytic in dehydrogenating adsorbed 
coronene than Cu(lll) in dehydrogenating benzene [15] . 
We therefore propose to use coronene as the first-step 
carbon source to seed carbon clusters on the patterned 
surface. 

After the deposition and dehydrogenation of coronene 
on Mn-Cu(lll) surface, all the adsorbed 7CR islands will 
have the same orientation. In particular, when two such 
7CR islands coalesce to form a larger graphene cluster, 
no GB is formed (Fig. 4). In order to achieve a continu- 
ous sheet of graphene, we invoke a second step of growth 
to supply C atoms to fill the openings between the car- 
bon clusters. We notice that, in contrast to Cu(lll), 
on which C adatoms are energetically much more favor- 
able to nucleate than to stay apart [16] . here on Mn- 
Cu(lll) the carbon adatoms are less strongly inclined 
to nucleate. The energy difference between a C dimer 
and two C monomers on the substrate, calculated as 

— -Emonomer+substrate X 2 -^substrate -^dimcr+substrate 

is 2.92 eV on Cu(lll) and becomes 1.70 eV on Mn- 
Cu(lll). Therefore, the conventional CVD growth us- 
ing methane or ethylene could be applied here to sup- 
ply carbon adatoms to diffuse and attach to the nearby 
coronene-seeded and correctly-oriented carbon islands 
(Fig. 4), rather than to nucleate new islands, which is 
similar to the enhanced layer-by-layer growth of Ag on 
Ag(lll) via a two-step kinetic pathway [25] ■ Eventually, 
the 7CR-seeded islands will be enlarged and connected to 
achieve a single-crystal graphene sheet with no or greatly 
suppressed GBs. 

It is important to note that, when a 7CR seed grows 
larger, the carbon clusters will still be at or close to the 
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FIG. 4: (color online). Enlargement of coronene-seeded car- 
bon islands via conventional CVD growth. The individual C 
adatoms supplied in the second step of the "seed and grow" 
kinetic pathway diffuse and attach to the nearby islands to fill 
the opening spaces, resulting in a larger graphene sheet with 
no GBs. 



HSO, because an edge C atom prefers to reside at the 
bridge site between a Mn and a Cu atom. Because all the 
islands have nearly the same orientation and dome-like 
geometry, they will be able to make minimal local ad- 
justments when they meet, and coalesce to form a single 
larger graphene sheet without GBs. Because of the lat- 
tice mismatch and its stronger interaction with the Mn- 
Cu(lll) surface than with pure Cu(lll), graphene may 
have a corrugated geometry, similar to that on Ru [41] . 
Finally, the playground of using patterned substrates 
is not necessarily limited to the Mn-Cu(lll) surface; 
other superstructured surface alloys with different tran- 
sition metals beyond the ones already considered here 
are also worth exploring. The present study of graphene 
growth on patterned substrates via a two-step kinetic 
process thus opens the door towards a new and viable 
approach for mass production of single crystalline mono- 
layer graphene. 
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